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Abstract 

The dynamics of charge qubits coupled to a nanomechanical resonator under influence of both 
a phonon bath in contact with the resonator and irreversible decay of the qubits is considered. 
The focus of our analysis is devoted to multi partite entanglement and the effects arising from 
the coupling to the reservoir. It is shown that despite losses, entanglement formation may still 
persist for relatively long times and it is especially robust against temperature dependence of the 
reservoir. Together with control of system parameters, the system may therefore be especially 
suited for quantum information processing. Furthermore, our results shed light on the evolution of 
open quantum many-body systems. For instance, due to intrinsic qubit-qubit couplings our model 
is related to a driven XY spin model. 

PACS numbers: 74.81.Fa, 03.65.Ud; 03.65.Yz. 
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I. INTRODUCTION 



The interface between quantum and classical mechanics has been an issue of intense 
research ever since the very early days of quantum physics. Until recently, the topic was 
restricted to the theoretical community, but the experimental progress seen during the last 
decade has made it possible to now study quantum properties of objects with a large number 
of degrees of freedom. For instance, the double slit experiment with C — 60 molecules has 
been realized^, and the decoherence of superposition electromagnetic coherent states with 10 
photons has been measured^ (which later was improved to a situation with 36 photons^). A 
system that allows for experimental studies of the classical limit is an optomechanical cavity. 
Here a cavity mirror or a membrane interacts via radiation pressure with the cavity field. At 
the moment, utilizing the interaction of the cavity field, the mechanical mirror /membrane 
can be cooled down to temperatures corresponding to a state very close to the ground 
state^. Yet another but related system is nanomechanical resonators, which has been cooled 
down to the mK regime by means of coupling to single-electron transistors. At such low 
temperatures, quantum properties crucially begin to play an important role^. Both these 
mechanical systems pave the way to the study of quantum phenomena, such as entanglement 
and the superposition principle, on a macroscopic level^. 

Apart from exploring the borders between classical and quantum mechanics through en- 
tanglement of macroscopic objects^, nanomechanical resonators provide also a way for multi- 
particle entanglement generation^, a subject that has attracted a great deal of interest lately, 
see^ and references therein. In particular, via piezoelectric interactions, nanomechanical os- 
cillators can be coherently coupled to many Joscphson junctions (JJs) simultaneously-'^^'^''^. 
Indeed, an important problem associated with quantum computation and information is 
the problem of engineering entanglement in multi-particle systems. At the same time as 
controllable multi-particle entanglement schemes are necessary for quantum information 
processing, equally important is the understanding how decoherence affects these systems. 
Lately it has been shown that interesting aspects may arise when a system is coupled to an 
environment, different from the general idea that the reservoir will solely induce decoher- 
ence and thereby destroy entanglement. For example, Kim et al. demonstrated that initially 
uncorrelated quibits coupled to a thermal field can become entangled^^. Later, Eberly and 
co-workers showed that the decay of entanglement shared between two initially entangled 
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qubits, due to coupling to a reservoir, can become sudden, so called entanglement sudden 
death^^. As entanglement might be lost instantaneously, it might also be regained suddenly, 
the entanglement sudden birth^^. 

Motivated by, at the one hand the experimental progress in recent years, and at the 
other because of the need for a deeper understanding of the multi-particle entanglement, 
we study the dynamics of a nanomechanical oscillator coherently coupled to charge JJs. 
Losses of both the quibits as well as the resonator are taken into account by solving the 
full system-reservoir master equation. We find both entanglement sudden death and birth 
dynamics, i.e. environmentally generated entanglement. Furthermore, our master equation 
approach allows for the investigation of temperature dependence, and we especially find that 
the sudden birth of entanglement is relatively insensitive to the reservoir temperature. This 
is further concluded from the fact that the many subsystems (JJs and resonator) disentangle 
at specific times regardless of the initial state of the resonator. Due to the effective dipole- 
dipole couphng appearing between the JJs, the model has resemblances with spin XX and 
XY models^^ and thereby gives insight into decoherent evolution of generalized interacting 
spin systems. We, however, restrict our analysis, due to exponentially growth of the Hilbert 
space dimensions, to a maximum of three JJs coupled to the nanoresonator. 

This paper is organized as follows. In Sec. II we derive the effective Hamiltonian of 
the system for arbitrary number of JJs. The general master equation is then introduced. 
Section III presents our results, where we in particular study concurrence and the tangle. 
The paper is concluded in Sec. IV. 

(Color online) Schematic diagram of N charge qubits coupled through the nanomechanical 
resonator. Here, each superconducting-quantum-interference-device loop (SQUID) connect- 
ing the superconducting island is represented by an effective JJ of coupling energy Ej and 
capacitance Cj. The Cooper-pair box is biased by a voltage Vg through the gate capaci- 
tance Cg and coupled to the bulk superconductors by two identical small JJs and 0c is the 
field-induced flux through the SQUID loop. 

II. MODEL SYSTEM 

To begin with, we consider a system of N charge qubits coupled to a nanoresonator via 
electrostatic interaction^^. We denote by and a the creation and annihilation operators 
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for a resonator phonon. Assuming weak or moderate vibration amplitudes n = a^a, and an 
energy splitting of the two lowest qubit states \g) and |e) to be resonant or quasiresonant 
with the single phonon energy, it is legitimate to restrict the analysis to a qubit sub-space 
only spanned by these two states^''. Such a reduction of the Hilbert space brings out an 
effective dipole-dipole coupling between the qubits^^. Turning to an interaction picture, the 
A'^-qubit Hamiltonian becomes^^ 

N 

j,i=l,i^j 

+ E l^^i'^ - l^J^^'^ +^{a^ + a) a^^]. (1) 
j=i ^ J 

Here, x is the effective dipole-dipole coupling strength, V the gate voltage, Ej is the coupling 
energy of the Josephson junction and 00 = 2(c'^+Cj) \/'^/2^C^, is the effective resonator-qubit 
coupling, where e is the electron charge, Ci;' is the capacitance parameter, which depends on 
the thickness of the junction, Cj is the junction capacitance, Cg is the gate capacitance which 
is externally controlled and used to induce offset charges on the island, v is the resonator 
frequency (see Fig. 1). The cr-operators are the standard Pauli matrices in the qubit basis 
\g) and |e); [o"q,,(3"/3] = 2i5Q,^^o"^, az\g) = —\g) and az\e) = \e). Note that we have further 
assumed that all parameters are independent on charge qubits, e.g. the qubits are taken 
identical and the gate voltage the same for each of them. 

In the absence of the nanoresonator, the above Hamiltonian is a special case of a driven 
XY spin model; the gate voltage has the meaning of an external field, x the spin-spin cou- 
pling and Ej a resonant driving of each spin^^. Furthermore, once coupled to the resonator, 
the model is related to the Dicke model of quantum optics where it has been shown that 
the dipole-dipole interaction induces novel phases such as spin glasses^^. 

To take resonator as well as qubit losses into account, we follow the standard approach 
by considering a Markovian master equation^°. This gives an equation of motion for the full 
system density matrix p, and takes the form^^ 

dp x 

— = - j^Iin, p] - {a) ap - 2apa) + pa^ a) 



1 ^ 
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E 7..^^ i^^^^P + P^^ - ^^?P^ ■ (2) 

The first term gives the unitary time evolution, while the remaining ones correspond 
to the effective couphng to the reservoirs. We symbohze by N, the number of thermal 
phonons at the resonator frequency v. Here we consider the case of a thermal bath without 
squeezing. In Eq. (2), k is the resonator decay rate, and the 7jj's are effective decay rates 
of qubit j, which will all be taken equal = 7. The qubits undergo as well a correlated 
decay with rates F = ^ j)i leading to further decoherence. 

In the representation p{t) — exp {\Ht) p{t) exp (— the equations of motion for the 
density matrix attain a simple form. For arbitrary initial states and using Eqs. (1) and 
(2), we numerically obtain the matrix elements of the density operator. For the special 
case of a system consisting of a single lossless qubit coupled to the resonator, one can 
obtain exact analytical forms of the diagonal matrix elements by imposing the rotating wave 
approximation in which rapidly oscillating terms are neglected. 

III. DYNAMICS 

We proceed by obtaining the solutions of the master equation introduced in the previous 
Section. One of the important aspects, especially for open systems, is the development of 
quantum correlations. The rule of thumb is that such correlations are quickly diminish in 
the presence of decoherence, but we will demonstrate that this is not necessarily the case. 
To this end, in this Section we present a detailed study of entanglement formation for two 
qubits coupled to a nanoresonator. 

For two-qubit systems, the concurrence is a measure of the entanglement for both pure 
and mixed states^^'^^. For a density matrix which represents the state of a bipartite 
system, the concurrence is defined as C{p) = max[0, Ai — A2 — A3 — A4], where the Aj's are 
the non-negative eigenvalues, in decreasing order (Ai > A2 > A3 > A4), of the Hermitian 
matrix T = -\/v^^V^ and p — {ay ay) p* {ay (8) ay). Here, p* is the complex conjugate 
of the density matrix p when it is expressed in a fixed basis and ay represents the Pauli 
matrix in the same basis. The function C{p) ranges from for a separable state to 1 for a 
maximum entanglement. The tangle in a pure two qubit system is defined as the square of the 
concurrence^^'^®. Rungta et al.^'^ defined the /-concurrence in terms of a universal-inverter 
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FIG. 1: (Color online) Time evolution of the tangle (3) between a single charge qubit and the 
resonator as a function of the scaled time cot and resonator losses k/uj. The J, J was initially in its 
excited state and the resonator in the vacuum state. The parameters are, Ej = V/oj = 1, and 
7 = r = 0.0. 

which is a generahzation to higher dimensions of the two spin flip operation. However, 
for arbitrary bipartite mixed states one should minimize the average squared pure state 
concurrence^*. In this case, the generalized mixed state form of the tangle can be written as 

r(p) = 2min5^(l-Tr[(p»n), (3) 

where pi*^ is the marginal state for the i^^ term in the ensemble decomposition and Si is 
a convex combination of the pure states. It is interesting mentioning here that, when we 
consider the two-qubit case, Eq. (3) gives the same results as the concurrence. 

In the first scheme, we assume that a single charge qubit, i.e. = 1, is coupled to the 
resonator which is damped at a rate while the qubit losses, on the other hand, is assumed 
zero 7 = r = and Nth= 0. In this case, we consider the qubit initially in its excited state 
|e) and the initial state of the resonator is |0). The evolved global state p(t) is determined 
after finding the time evolution of each matrix element, by using Eq. (2). The resonator and 
qubit are two interacting systems, since the ranks of their reduced density operators {pr,q(t)) 
are not greater than two, Pr,q(t) = Tr^ ,. {pit)}. Hence, the use of tangle (concurrence) as a 
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FIG. 2: (Color online) Time evolution of the tangle (3) between two charge qubits as a function 
of the scaled time cot and resonator losses k/uj. The first qubit is initially in its excited state, 
the second qubit in its ground state, and the resonator in the vacuum state. The parameters are, 
Ej = u, V/lo = 1, X = 0, and 7 = r = 0.0. 
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FIG. 3: (Color online) Tangle between two charge qubits as a function of the scaled time cot. The 
initial state of the resonator is a mixed state (see text for specific parameters) with different values 
of the decay rate 'y/co, where 7/0; = 0.01 (solid red curve), 7/0; = 0.1 (dotted green curve) and 
7/0; = 0.7 (dashed blue curve). The other parameters are x = 0; ^ = 0; ^th = 0.5, P/w = 0.001, 
and Ej = u/uj. 
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FIG. 4: Tangle between two charge qubits as a function of the scaled time ujt. The initial state 
is IV'(O)) = |ei,e2;0) with different values of x/^^ where (a) x/^ = 30, (b) x/^ = 15 and (c) 
x/uj = 0.01. The other parameters are the same as Fig. 3, but here, k = and 7 = F = 0.0. 

measure of quantum correlations/entanglement between the two subsystems applies in this 
case. In order to calculate the tangle, we choose the standard basis of the product space 
of each subsystem as |0,e), \0,g), |l,e) and (here "1" indicates the one phonon Fock 

state), then one can effectively consider the density matrix of the system formed by two 
two-level subsystems at each instant of time t. 

Figure I displays the tangle with dynamics especially showing a sequence of entanglement 
sudden deaths/births. Its evolution is strongly affected by the modulation term that depends 
on the decay rate n. As the resonator decay rate increases, this effect is less pronounced 
and as expected, in the bad resonator limit all entanglement gets rapidly lost. It is seen 
that in the k, —>■ limit we still have entanglement sudden deaths/births lasting over a 
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large time span. Such features might seem remarkable since the von Neumann entropy (not 
presented here) in the rotating wave approximation is indeed a smooth function of time t 
in the k, = limit. It is known that even without the rotating wave approximation, the 
von Neumann entropy has a smooth time evolution^^. To clarify this point, we stress the 
fact that, the von Neumann entropy measures entanglement of the reduced density matrix, 
while the tangle r (p), used in this work, utilizes the evolved global state p{t) together with 
a reduced phonon basis in which the system substitutes two two-level subsystems. These 
aspects imply that the von Neumann entropy and the tangle measure different properties 
and thereby give different behaviors; one being smooth and one possessing discontinuous 
time derivatives. 

In general, the time evolution is quite complicated. However, within the rotating wave 
approximation and the special case of x = 7 = F = 0, taking into account the initial state 
of the system is p(0) = |e, 0)(0, e|, the tangle takes a simple form. 



T (p) = max 



0,sin(2a;i)|2cosh^— + cosh^— c^^t^ 



— 2sinh (^—cut^ — sinh (—out 



(4) 



In getting Eq. ( 4), we consider the temperature T = K, i.e. the average thermal phonon 
number is zero. This expression verifies the fact that entanglement vanishes in the large 
decay limit k — > 00. 

We now turn to the case of two qubits, opening up for multi partite entanglement. 
Thereby, since the qubits possess an indirect dipole-dipole interaction among each other, 
there is an additional term that drives the entanglement evolution. In the two qubit case, 
we assume the initial state of the qubits as 



p"(0) = a|e,e)(ee| +6|e,^)(e^| + (1 -a)|^^)(^c/| 

+f\e,g){ge\+r\ge){eg\+c\ge){ge\, (5) 

where a, b, c and / are constants. To compare with the single qubit case, we assume a — 1 
and b — c — f — 0, i.e. both qubits start out in their excited states. Prom Fig. 2, we 
see that the problem of entanglement creation and evolution in the two-qubit system is not 
similar to the one-qubit case. In the present case, the entanglement does not show a zero 
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FIG. 5: (Color online) Tangle as a function of the scaled time Lot for the case with two initially 
excited charge qubits and different values of N, where N = 0.01,0.1 and 0.5, (dashed red, dotted 
green, and solid black lines respectively). The other parameters are the same as Fig. 3 and k = 0. 

value lasting for long times (compare Fig. 2 and Fig. 3). Also, for small values of the decay 
rate k, additional drops in the tangle are seen. However, the entanglement vanishes in the 
large decay limit in a similar way as for the single qubit case. 
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FIG. 6: Tangle as a function of the scaled time Lot, (for the three qubit case). The initial state of 
the resonator is a thermal state with N = 0.5,. The other parameters are the same as in Fig. 6. 
The initial qubit state is taken as \eee). 



The time evolution of the tangle between two charge qubits for various values of the 
effective decay rate of the qubit (7/0;) is displayed in Fig. 4. Here it is assumed that the 
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initial state of the qubits is a mixed state, i.e. a — 0.1 and 6 = c = 1 = 1. In this figure, 
we see that the entanglement is sensitive for 7/c<j, but for sufficiently small values of 7/0; it 
is relatively robust and persists over longer time periods. However, for large values of 7/0;, 
the sudden death of entanglement occurs. 

In Fig. 5, the time evolution of the tangle between two qubits is plotted against the scaled 
time cut for nonzero dipole-dipole coupling strengths (x 7^ 0) when the qubits are initially 
in the state 1 61,62). In this figure we ignored the effect of losses k, and 7. It is observed 
that the entanglement creation time decreases as x becomes large. Which means that, this 
parameter can be used as a controller for entanglement creation, in agreement with You, 
et aP° where they stated that this type of spin-spin coupling can be conveniently used to 
formulate an efficient quantum-computing scheme. 

When the mean thermal phonon number is not zero, entanglement shows completely 
different time behavior. This is illustrated in Fig. 6, where we plot the tangle as a function 
of the scaled time for different values of the mean thermal phonon number. In contrast to 
the previous cases where entanglement was created immediately after t = 0, we note that 
the creation of entanglement is now delayed, that it occurs after a finite time. The reason 
for the delayed creation of entanglement can be understood as follows. When the system 
is prepared in the state |6i,62), the resulting transitions are cascades: The system decays 
first to the intermediate states, from which it then decays to its ground state. Since the 
transition rates to and from the intermediate states are different when N ^ 0, there appears 
unbalanced population distribution between these states. This may result in a transient 
entanglement between the qubits. Also, we see that in this case the entanglement sudden 
death always happens no matter which entangled state the qubits are initially in and no 
matter how small the nonzero mean thermal phonon number is. This implies that, even for 
relatively "hot" baths the system shows entanglement properties. Similar results were found 
in a system of two lossless two-level atom interacting with a finite Q cavity^^. 

Another important question one may ask: Docs the properties of entanglement sudden 
birth depend on the number of qubits {N > 2), and especially is the effect less pronounced for 
more qubits? To answer this question, we plot in Fig. 7 the tangle with respect to the scaled 
time for three qubits (A^ = 3). This figure shows over a range of interaction times that the 
entanglement start to build up at earlier times closer to the initial interaction time. This is 
in striking contrast to the two qubits case displayed in Fig. 6 where the entanglement starts 
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to build up at a much later time. Also, the number of sudden deaths/births is increased 
with less amphtudes. This might be interpreted as the state's entanglement becoming more 
robust when the size of the system increases. However, what really matters is not that the 
initial entanglement does not disappear fully but that a significant fraction of it remains, 
either to be directly used, or to be distilled without an excessively large overhead in resources. 
Overall, we find that the number of qubits acts as well as a control parameter for the qubit- 
qubit entanglement. 

IV. CONCLUSIONS 

In this paper, we have investigated the time evolution of the entanglement between a 
nanomechanical resonator and N — 1, 2, 3 charge qubits. We have shown that by changing 
the dipole-dipole coupling, which can be achieved by varying the mutual distances between 
the JJs, one can control the entanglement. The decay which is regarded as a source of 
decoherence, creates in general entanglement sudden birth followed by a sudden death. Thus, 
entanglement may in fact be generated via the environment. Finally, the entanglement 
of three qubits is discussed. We have in particular demonstrated that by increasing the 
number of qubits, the entanglement starts at earlier times and shows more frequent sudden 
deaths/births. 

The resilient properties of the system together with its entanglement robustness therefore 
makes it a suitable candidate for quantum information processing, fundamental studies of 
quantum many-body dynamics, or analyzes of the quantum-classical interface. Its nanome- 
chanical structure might as well have applications in precision measurements^^. 
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